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Abstract 

We review the recent progress on the numerical implementation of the Loop-Tree Duality Method (LTDM) for the 
calculation of scattering amplitudes. A central point is the analysis of the singularities of the integrand. In the frame¬ 
work of the LTDM some of these singularities cancel out. The ones left over are dealt with by contour deformation. 
We present details on how to achieve this as well as first results. 
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1. Introduction 

When calculating NLO (NNLO) cross-sections one 
needs to consider the tree- and loop-contributions sepa¬ 
rately. Especially loops with many external legs prove 
to be challenging. Considerable progress has already 
been made in order to attack this problem; OPP- 
Method, Unitarity Methods, Mellin-Barnes Representa¬ 
tion, Sector Decomposition III. The advantage of these 
methods is that they made possible what was impossible 
before, but still a lot of effort has to be put in to cancel 
infrared singularities among real and virtual corrections. 
Additional difficulties arise from threshold singularities 
that lead to numerical instabilities. The Loop-Tree Du¬ 
ality method aims towards a combined treatment of tree- 
and loop- contributions. Therefore the Loop-Tree Dual¬ 
ity method casts the virtual corrections in a form that 
closely resembles the real ones. 


2. Loop-Tree Duality at one loop 


The most general, dimensionally regularized one-loop 
scalar integral can be written as 121 : 
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with the Feynman propagator Gfiqf) - [qf - m? + , 

internal momenta qt - (\ + p\ + ■ ■ ■ + pi = {\ + k, and 
shorthand integral notation = -i Jd‘^{il{ 27 :Y . As 
a first step, one performs the integration over the com¬ 
plex energy components of the loop four-momentum by 
applying the residue theorem. The integration contour 
is chosen such that it encloses the poles with positive 
energy and negative imaginary part, see figure below; 



Figure 1: Location of poles and integration contour Cl in the complex 
f’l o-plane. 


The residue theorem is employed by taking the 
residues of the poles inside of the contour and summing 
over them. Given an appropriate gauge choice the in¬ 
tegrand in eq. Q contains only simple poles. Thus 
the residue of an individual pole is done by taking the 
residue of a single propagator and evaluating the other 
propagators at the position of the residue 
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= {\GD{qi\qj). (3) 
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The subscript “+” on the right hand side of eq. 0 in¬ 
dicates that the positive-energy solution is to be taken. 
Furthermore 77 is a future-like vector, i.e. rf >Q,rio> 0. 
It is dependent on the choice of coordinate system, how¬ 
ever it cancels out once one adds all dual contributions. 
Hence physical objects like scattering cross sections 
will stay frame-indepedent. Evaluating the “non-cut” 
propagators at the position of the pole leads to a mod¬ 
ification of the usual Feynman prescription. In eq. 0 
it is shown that instead one ends up with the so called 
“dual prescription” which serves to keep track of the 
correct sign of the iO-prescription of the corresponding 
propagator. Collecting all the pieces and putting them 
together, one arrives at 


• The singularities of the loop diagram appear as sin¬ 
gularities of the Dual Integrals. 

• Tensor loop integrals and physical scattering am¬ 
plitudes are treated in the same way since the 
Loop-Tree Duality works only on propagators. 

• Virtual corrections are recast in a form, that closely 
parallels the contribution of real corrections. 

This is the formalism for the one-loop case. Solutions 
for more complicated situations like multiple loops 13 
or higher order poles @1 are described in the respective 
references. 

3. Singular behavior of the loop integrand 

As a preparatory step it will prove useful to introduce 
an alternative way of denoting the dual propagator. This 
will give a more natural access to its singularities. 


r - ^ 

L^^\PuP2,---,Pn)^-'^ I j~[ Gofe; ?;•) 

7=1 

j*i 

(4) 

with 6(qi) - 2m6+(qJ - mf). Thus, by virtue of em¬ 
ploying the residue theorem, it is possible to rewrite a 
one-loop amplitude as a sum of single-cut phase-space 
integrals over the loop-three-momentum. The i-th dual 
contribution has the i-th propagator set on-shell while 
the left over Feynman propagators get promoted to Dual 
propagators. 

Lf^\pi,P2,---,PN) = 

Pi Pi+l 



MqdGoiqhqj) = 2.711 
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with kji - qj - qi and qG^ - -I- /«? - lO. 

In fig. the on-shell hyperboloids of three propaga¬ 
tors in loop-mometum-space are sketched. 



Figure 2: Graphical representation of the solution of the LTDM at 
one-loop. 


Figure 3: On-shell hyperboloids for three ai'bitrary propagators in 
Cartesian coordinates. 


The LTDM features a couple of interesting proper¬ 
ties: 

• Number of single cut Dual Contributions equals 
the number of legs, this way a loop diagram is fully 
opened to tree diagrams. 


The loop integrand becomes singular at hyperboloids 
with qG^ - + inf - iO (solid lines) and qG^ - 

qi^ + inf - iO (dashed lines) and origin in Ap¬ 
plying the LTDM is equivalent to integrating along the 
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forward hyperboloids. The intersection of two for¬ 
ward hyperboloids (solid lines) leads to singularities 
that will cancel among dual integrals. The intersection 
of a forward with a backward hyperboloid (solid line 
with dashed line) leads to a singularity that remains and 
therefore has to be dealt with by contour deformation. 
Singularities appear where the denominator of the dual 
propagator becomes 0. Since the denominator in eq. 0 
has been rewritten as a difference of squares, one can 
use the third binomic formula to take it apart and ex¬ 
tract the conditions for which the singularities show up: 



Eq. 0 describes an ellipsoid in the loop three- 
momentum and demands < 0. An ellipsoid is the 
result of the intersection of a forward with a backward 
hyperboloid. The origins of the hyperboloids are sepa¬ 
rated in a time-like fashion, expressed by the conditions 

kj,- - (nij + iriif > 0, kjifi < 0. (8) 

The singularity described by eq. 0 has a hyperboloid 
shape as a result of the intersection of two forward on- 
shell hyperboloids of space-like separation. The corre¬ 
sponding condition is 

- ('«./■ - ^ 0 - ( 9 ) 

Here, ky, o may be positive or negative. 

4. Numerical Implementation 

The singularities associated with ellipsoid intersec¬ 
tions require contour deformation. The following 
one-dimensional example illustrates how this can be 
achieved. The function 


exponent is a large negative number, hence exponenti¬ 
ating it suppresses the deformation. 


Shape of the contour deformation 



Figure 4: Contour deformation as in eq. GI} for £ = 2 and A = 5. 


For the integration in loop three-momentum-space 
eq. ([nj needs to be generalized to three dimensions. 
This is done by modifying the exponent and, of course, 
promoting to 


/ ^ / H- /d/exp- — 

width 
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GJ' now plays the same role as did in eq. 0- 

The parameter “width” allows to control the width of 
the deformation, thus giving additional control over it. 



Figure 5: Hyperboloid and ellipsoid singularities in a triangle. 


/( 4 ) = 


1 

£l-E^ + iO 
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has poles at ±(E - iO). Simply integrating along the real 
axis would lead to infinities. Therefore the integration 
contour has to be deformed to go around the poles. A 
suitable contour deformation would be: 

, / el-E^\ 


The parameter A serves to scale the deformation along 
the imaginary axis. At the position of the pole the expo¬ 
nent becomes 0 and thus the exponential function hits 
its maximum, which is 1. Far away from the poles, the 


The example in fig. 0 shows the location of the sin¬ 
gularities as they might appear for a triangle. Every box 
corresponds to a Dual contribution, i.e. cut. The easiest 
case is the third box, which doesn’t display any singu¬ 
larities of any type and thus can be integrated straight¬ 
forwardly. In the first contribution there is only a hyper¬ 
boloid singularity, which cancels with the hyperboloid 
singularity from the second contribution. But since the 
second contribution features an ellipsoid singularity as 
well, this contribution requires contour deformation. In 
turn, in order to preserve the cancellation of the hyper¬ 
boloid singularities, contribution one will need the ex¬ 
act same deformation, despite no ellipsoid being present 
there. This means that hyperboloid singularities act as a 
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coupling between different dual contributions. The cou¬ 
pled contributions need to be deformed according to all 
ellipsoid singularities which they both have. 

In the massless case collinear singularities arising from 
forward-forward intersections cancel among Dual in¬ 
tegrals, because the on-shell hyperboloids of fig. 0 
would be tangential there. Nonetheless, collinear and 
soft singularities from forward-backward intersections 
remain, but they are restricted to a finite region that can 
be mapped to the real phase-space emission Q. 

5. Results and Conclusion 

A computer program has been written in C-n-. The nu¬ 
merical integration step is performed with Cuhre from 
the Cuba-Library 0. Analytic values for comparison 
are generated with LoopTools (LT) Q. The code runs 
on an Intel i7 desktop machine with 3.4GHz. It is capa¬ 
ble of calculating triangles, boxes and pentagons with 
no deformation needed with 4 digits precision in 0.5s. 
The table below shows an explicit result where a pen¬ 
tagon with deformation has been calculated with 4 dig¬ 
its precision in 25s time. 



Real part 

Re Err. 

Imaginary part 

Im Err. 

LT 

-1.001066E-10 

0 

-5.208136E-10 

0 

LTDM 

-1.001089E-10 

9.1E-16 

-5.208556E-10 

9.1E-16 


Another test the program has been put through is a 
scan of the region around threshold. To produce the two 
graphs in fig (|^ and (|^, a triangle has been taken; the 
center-of-mass energy s has been kept constant while 
the mass m (which was taken to be identical for all three 
internal lines) was varied. 



Figure 6: Real part. The curve is LoopTools, the dots are the LTDM. 



Figure 7: Imaginary part. The curve is LoopTools, the dots are the 
LTDM. 


From the plots one can see that LoopTools’ analytical 
values are well matched by the computer program, this 
indicates good precision. There is no loss of precision 
or increase in calculation time close to threshold. 
Despite still being in development, the program already 
shows competitive speed and precision when calculat¬ 
ing loop integrals. At the moment, only diagrams up to 
pentagons have been checked, but within the framework 
of the LTDM, extension to graphs with more external 
legs is straightforward and easy to realize. 
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